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ABSTRACT 

We formulate a general path integral approach which describes statistics of current fluctuations in meso- 
scopic coherent conductors at arbitrary frequencies and in the presence of interactions. Applying this 
approach to the non-interacting case, we analyze the frequency dispersion of the third cumulant of the 
current operator £3 at frequencies well below both the inverse charge relaxation time and the inverse elec- 
tron dwell time. This dispersion turns out to be important in the frequency range comparable to applied 
voltages. For comparatively transparent conductors it may lead to the sign change of S3 . We also analyze 
the behavior of the second cumulant of the current operator <So (current noise) in the presence of electron- 
electron interactions. In a wide range of parameters we obtain explicit universal dependencies of 1S2 on 
temperature, voltage and frequency. We demonstrate that Coulomb interaction decreases the Nyquist noise. 
In this case the interaction correction to the noise spectrum is governed by the combination Y2 n T n (T n — 1), 
where T n is the transmission of the n-th conducting mode. The effect of electron-electron interactions 
on the shot noise is more complicated. At sufficiently large voltages we recover two different interaction 
corrections entering with opposite signs. The net result is proportional to Y] T n {T n — 1)(1 — 2T„), i.e. 
Coulomb interaction decreases the shot noise at low transmissions and increases it at high transmissions. 
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1. INTRODUCTION 

Experimental and theoretical studies of mesoscopic conductors reveal a rich variety of low temperature 
properties and effects caused by interplay between scattering, interactions, quantum coherence and electron 
charge discreteness. Investigations of the whole scope of these effects are of primary interest because of 
their fundamental importance as well due to rapidly growing number of their potential applications. 

A great deal of information is usually obtained from studying electron transport. Additional information 
can be extracted from investigations of fluctuation effects. During last years much attention has been 
devoted to the shot noise 1 described by the second moment of the current operator. One can also study 
higher order correlators of the current operator thereby extending the amount of information already 
obtained from electron transport and shot noise. Recently the first experimental study of the third current 
cumulant in mesoscopic tunnel junctions was reported. 2 

A theoretical framework which enables one to analyze statistics of charge transfer in mesoscopic con- 
ductors was developed in Ref. 3. This theory of full counting statistics (FCS) allows to evaluate any 
cumulant of the current operator in the absence of interactions and in the zero frequency limit. Under 
these conditions higher order current cumulants were investigated by a number of authors. 4-7 In order to 
include interactions and to analyze frequency dispersion of current fluctuations it is necessary to go beyond 
the FCS theory and to develop a more general real time path integral technique. 8-12 

The goal of the present paper is to address statistics of current fluctuations at non-zero frequencies and 
in the presence of electron-electron interactions. In what follows we will develop a general path integral 
approach which allows to perform a complete analysis of electron-electron interaction effects in mesoscopic 
coherent conductors described by an arbitrary - though energy independent - scattering matrix. This 



approach also provides a straightforward generalization of the FCS theory 3 to non-zero frequencies. Of 
central importance in this context is a general and formally exact expression for the real time effective 
action of a coherent conductor. This expression will be derived and analyzed in various limiting cases in 
Section II of the paper. With the aid of our technique in Section III we will investigate the frequency 
dispersion of the third cumulant of the current operator in mesoscopic coherent conductors. Section IV is 
devoted to the analysis of the electron-electron interaction correction to the second cumulant of the current 
operator, i.e. to the current noise. The paper is briefly concluded in Section V. 



2. THE MODEL AND EFFECTIVE ACTION 
2.1. General analysis 

In our analysis we will use the real time path integral formalism developed for the systems of interacting 
fermions. After the standard Hubbard-Stratonovich decoupling of the interaction term in the Hamil- 
tonian one can exactly integrate out fermions and arrive at the effective action S which depends on the 
fluctuating fields Vi^ 2 (t,r). Let us define 
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with the trace taken over the fcrmionic variables. Here po is the initial N— particle density matrix of 
electrons, 
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Here Hi 2 are the effective Hamiltonians on the forward and backward parts of the Keldysh contour, U (r) 
describes the static potential, T and T are respectively the forward and backward time ordering operators. 
Integrating out fermions in Eq. (1), we obtain 



where 



iS = 2Trln[l + (u 2 ux — l)po], 
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are the evolution operators pertaining to the Hamiltonians (2), po stands for the initial single-particle 
density matrix. 

The total action of the problem S also includes the contribution of the fluctuating electromagnetic field 
outside the scatterer. This contribution reads 
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Here Zs(lo) is the effective linear impedance of the electromagnetic environment, and ips is the quantum 
phase variable related to the fluctuating voltage Vs across the the external leads by means of the equation 
ips = eVs. The superscripts (±) just reflect the fact that the phase variable is defined on the Keldysh 
contour, so that symmetric (^5) and antisymmetric (fg) combinations of the phases on two branches of 
this contour should be introduced. 

Although the above expression is valid for arbitrary Zs(uj), for the sake of definiteness and simplicity 
below we will choose 1/Zs(u>) = —iuuC+ l/Rs, where C is the barrier capacitance and Rs is the resistance 
of external leads (see Fig. 1). The total action of our system then reads 



S = S C 



S . 




Figure 1. Schematics of the system under consideration. The scatterer is placed in-between two big metallic 
reservoirs and is characterized by the Landauer conductance 1/7? and the effective capacitance C. The reservoirs 
are connected to the voltage source via leads with the total Ohmic resistance Rs- 

In order to evaluate the evolution operators u\ t 2 and derive the expression for the action (3) it is necessary 
to specify the model of a mesoscopic conductor. Here we will adopt the standard model of a (comparatively 
short) coherent conductor placed in-between two bulk metallic reservoirs, see Fig. 1. The electron dwell 
time td is supposed to be shorter than any relevant time scale in our problem. Energy and phase relaxation 
times are, on the contrary, assumed to be long, i.e. inelastic relaxation is allowed in the reservoirs but not 
inside the conductor. Under these assumptions electron transport through the conductor can be described 
by the energy independent scattering matrix 
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and the standard Landauer formula for the conductance 1/R = (e 2 /ir) T n = (e 2 /7r)Trr t can be used, 
where T n defines the transmission for the n-th conducting mode. 

The effective action (3) can be expressed via the fluctuating phase fields tpi^ which are in turn related 
to the jumps of the fields Vifi across the scatterer as tfii^ = &{Via,i — Vm^), where Vl,r are fluctuating in 
time but constant in space fields in the left and right reservoirs. We note that in this case the right hand 
side of Eq. (1) differs from the FCS generating functional 14 only by a gauge transformation. 

Within the above model the evolution operators were derived in Refs. 8 and 9. Combining these 
expressions with Eq. (3) after some algebra we find 10 
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Here we introduced ip + = (tp\ + (f2)/2, ip — tpi — (f2 and 
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Taking the trace in Eq. (7) implies convolution with respect to internal time variables (x and y). The total 
time span is denoted by t. 



Eq. (7) defines a formally exact effective action for a coherent conductor described by an arbitrary 
energy independent scattering matrix (6). This expression allows to fully determine statistics of current 
fluctuations at arbitrary frequencies and in the presence of interactions. In the case of equilibrium fluc- 
tuations the formula (7) represents a real time analogue of the effective action 15-17 derived within the 
Matsubara technique. A formula similar to Eq. (7) was also presented recently in Refs. 11 and 12. In 
addition we note that, provided the field (p~ does not depend on time, Eq. (7) coincides with the generating 
function considered, e.g., in the problem of adiabatic pumping through mcsoscopic conductors. 18 

In order to proceed with further calculations it is useful to investigate the properties of the general 
expression (7) and - whenever it is possible - to find its simplified representations. This goal can be 
achieved with the aid of several different approximations which can be appropriate in the corresponding 
physical limits. Some of these approximations and the resulting expressions for the effective action are 
specified below. 

2.2. Expansion in T n 

Provided channel transmissions T n of the scatterer are sufficiently low it can be useful to employ a regular 
expansion of the action (7) in powers of T n . The first order terms of this expansion just yield the standard 
AES action for tunnel barriers. 16 Here we will proceed further with this expansion and establish the second 
order in T n contributions to Sq. In order to recover all such terms it is necessary to expand the logarithm 
in Eq. (7) up to the fourth order. After a straightforward - although somewhat tedious - calculation we 
obtain 
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Here we introduce the dimensionless conductance of the scatterer g = 2 T n as well as its Fano factor 
j3 = 5^ n T n (l — T n )/ "52i n T n . Eq. (9) represents a complete expression for the effective action valid up to 
the second order in the transmissions T„. This expression involves no further approximations and fully 
accounts for the non- linear dependence on the fluctuating phase fields y ± . 

2.3. Reflectionless barriers 

Another physically important limit is that of reflectionless barriers r = 0. In this limit the general expression 
(7) can be reduced to a very simple form by means of the exact procedure which we outline below. The 
derivation is based on the following property of the Fermi function n{E) — 1/(1 + cxp(E/T)) : 

n(E)n(E + w) = i (l - coth ^) n(E) + 1 (l + coth |L) n(E + u). (10) 



This relationship plays the same role as the Ward identity in the Luttinger liquid theory. Rewriting (10) 
in the coordinate representation one finds 

p (x - y)p {y -z) = (f(x - y) + f(y - z))p {x - z), (11) 

where f(x) = J ^ \ (coth ^ + l) c~ iuJX = -ifV coth[7rTx] + \8(x). Let us formally expand the action 
(7) in powers of the transmission matrix vt. This expansion gives rise to the integrals of the type 

J N = I dx l dx 2 - ■ - dx N a(x 1 )p (x 1 -X2)a(x 2 )po(x2-x 3 )a(x 3 )...pa(x N „ 1 -x N )a(x N )p (x N ~x 1 ), (12) 



where a(x) = e"^ (*) - 1. The identity (11) permits to evaluate such integrals exactly. It will be 
convenient for us to carry out N — 1 permutations of the variables Xj as follows. We start from the set (1) 
x\, x 2 , X3 . . . xn, and further introduce the following sets of integration parameters 

(2) xi x N x 2 x 3 ... x N -i 

(3) x\ x 2 x N x 3 ... x N -i 

(N-l) x\ x 2 ... x N - 2 x N Xn-1, 

i.e.. we successively insert xn in-bctween Xi and x 2 , x 2 and x 3 , an so on. Symmetrizing the integral (12) 
with respect to these TV — 1 permutations and excluding x^ from the functions po by virtue of the identity 
(11), we get 
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Eq. (14) allows one to relate Jm and Jn-i- After N — 2 iterations one finds the following result: 
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Now we are in a position to evaluate the trace A = Trln [l + 6(t — x)6(x)a(x)po(x — y)] . Exploiting Eq. 
(15) we obtain 
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Now we recall that in the case of refiectionless barriers (T n = 1 for all n) the action (7) reads 

iS = 27V ch {Trln [l + 9(t - x)6(x)a(x)p R (y,x)] + Trln [l + 8(t - x)6(x)a*(x)p L (y,x)] }, (17) 

where A^ c h is the total number of fully open channels in the conductor and Pl,r(v, x) = e l ^ L - R ^~ Vh - R ^ po{y- 
x). Expanding this action in powers of a(x) and a*(x) and applying the formula (16) we observe that the 

factors e PL < R ^~' f ' L > R ^ y ' cancel out exactly in all orders except in the first term of Eq. (16). This remaining 
term should be treated with sufficient care having in mind the formal divergence po(x,y) —> oo &t x ^ y. 
Finally we obtain 

iSo = -^- f dx^ + (x)W( V -{x)) + 9 - f dxdy , T l -W{<p-(x))W{<p-(y)), (18) 
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where g = 2iV c h and W(ip~) is the 2ir— periodic function equal to W(ip~) = ip~ for — it < Lp~ < it. Note 
that for the values of (p~ within the latter interval Eq. (18) exactly coincides with the action for an Ohmic 
conductor. Finally we point out that the existence of the exact solution (18) in the limit f = is not 
at all surprising since in that limit the problem can also be exactly treated by means of the bosonization 
technique. 19 



2.4. Local in time part of the action 

The above results for Sq demonstrate that the effective action can be split into two parts, 

So = S'local + •S'nonlocal, (19) 

determined respectively by local and non-local in time expressions. One can also observe that the term 
^locai is linear in ip + . Such splitting can be performed for arbitrary transmission values. It turns out that 
the local part of the effective action can be found exactly in all orders in T n . 

In order to derive this exact expression let us first fix ip + = eVt and tp~ = const. It is easy to 
demonstrate 10 that under these restrictions the action 5*0 reduces to the FCS generating function 3 
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where Arccosh a; = ln(a; + \J x 1 — 1) and v = eV/2T. The part S'local can be identified from this generating 
function in the limit DCl. Keeping only linear in v terms one can rewrite Eq. (20) in the form 
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Since Si oca i should be determined by the local in time combination one can now replace eV by an arbitrary 
function of time tp + and also allow for arbitrary changes of (p~ in time. This observation immediately 
brings us to the final result 
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Unfortunately the above exact expression for S'local alone is not sufficient for a full description of quantum 
properties of our system and no such simple derivation procedure exists for the non-local in time part 
S'nonlocal- Nonetheless, Eq. (22) provides an important piece of information about the properties of the 
effective action So and it can be successfully employed in further calculations. 

2.5. Expansion in ip~ 

In a large number of physical situations it is sufficient to restrict the analysis to the limit of large scatterer 
conductances j > 1. Either in this limit or, equally, in the case gs = 2ir/e 2 Rs 3> 1, fluctuations of the 
quantum field tp~ are considerably suppressed almost at all energies/frequencies. In this case it is useful 
to expand the exact expression for So in powers of </?~ keeping the full non-linearity in ip + in all orders 
of this expansion. Below we will present the result of this expansion up to the third order in tp~ . The 
corresponding derivation has been worked out in Ref. 9. It yields 
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where we have defined 7 = ^ T%{1 — T n )/J2 T n - Obviously, the local in time part of Eq. (23) can also be 
recovered by a direct expansion of Eq. (22) in powers of ip~ . 

Below we will make use of the result (23) in order to evaluate the frequency dependence of the third 
cumulant of the current operator and to derive the interaction correction to the shot noise spectrum. 

3. FREQUENCY DEPENDENCE OF THE THIRD CURRENT CUMULANT 

To begin with, let us recall that the standard way to describe the current noise in our system is to evaluate 
the symmetrized correlation function 

s(t,t') = ±{i(t)i(t') + i(t')i(t))-(if, (23) 

where / is the current operator in the circuit of Fig. 1. After the Fourier transformation within our path 
integral technique the same quantity can be expressed as follows 

*M = tre-J p ^ __g_ c^.-i (24) 

Here the path integral is taken over all possible configurations of the phases ( / 3 ± . In this way the effect of 
Coulomb interaction on current noise is fully accounted for by Eq. (24). This effect will be studied in the 
next section. 

Proceeding along the same lines one can also define higher cumulants of the current operator. Below we 
will only analyze the frequency dependence of the third cumulant. In doing so, we will ignore interactions. 
In this case the path integral is dominated by the trajectories <p + (t) = eVt and ip~ = 0, the latter condition 
being imposed after taking the functional derivatives: 
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It is important to emphasize that the current through the scattcrcr should be treated as a quantum operator. 
Since in general the current operators do not commute if taken at different times, the question arises about 
the proper ordering of such operators in the third and higher order correlation functions. One can verify 
that the definition (25) corresponds to the following combination of the correlation functions: 

S 3 (h,t 2 ,t 2 ) = l{(/(t x )(Ti(t 2 )/(*3))> + <(tJ(t 2 )/(t 3 ))% 
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- \(i(h))(f(ti)i(t 2 ) + i(t 2 )i(h)) + 2<%))<i(i 2 )>(/(t 3 )>. (26) 

This combination is of interest in the light of possible experimental investigations of current fluctuations. 
Indeed, the correlation function £2 is important because the symmetric combination of voltages V + can 
be viewed as a classical, measurable, voltage. 14 The noise is deduced from the measurable product 
V + (ti)V + (t 2 ), which is related to the symmetric correlator S 2 (23) or (24). Analogously, the measured 
product V + (ti)V + (t 2 )V + (£3) is related to the correlation function of the current operators S 3 defined in 
Eq. (26), see also Refs. 14 and 20 for further discussion of this point. 

Let us substitute the above expressions into Eq. (23) and, after taking derivatives over ip~ , set <p + (r) = 
eVr and ip~ — > 0. This is sufficient provided the time differences |£i — £ 2 |, |£i — £3! exceed the charge 
relaxation time trc and provided eV -C 1/trc". Eq. (23) then yields 
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where /(r) = Tsin(eVT/2)/sinh(7rTY). Performing the Fourier transformation 
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we arrive at the final result 
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Here we defined w = eV/T, Wi^ = ^i, 2 /2T and 

XH = cosh 2 cj + sinh 2 (u/4). (30) 

Eqs. (28)- (30) represent the main result of this section. They fully describe the third cumulant of the 
current operator at voltages and frequencies smaller than both 1/trc and 1/td. 

Let us briefly analyze Eqs. (28)-(30) in various limits. For u>i >2 = we recover the well known result 5 
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which in turn yields F — ► 1 in the limit of large voltages »> 1. Eq. (31) also holds for wi >2 <C v. 
In the limit v <C 1 one finds 

„. . 9 sinh w + sinhfSty) — 12iy cosh w , .„„. 

F(v<£l,w,0)=F(v<g.l,w,-w) = i — ( v 2 , (32) 

48 sinh w 

„, , , sinh(4w) + 4 sinh(2u0 — 8w cosh(2u») — 4w 9 

F(v < l,w, w) = — 5 — v. (33) 

128 sinh w cosh w 

These equations demonstrate that at large frequencies w 1 the function F decays exponentially with w. 
From Eqs. (32) and (33) we find respectively F oc v 2 e~ 2w /3 and F oc v 2 e~ 4w . 

Finally let us turn to the most interesting limit of low temperatures, in which case one always has 
v,Wi,W2 3> 1. Neglecting small corrections ~ l/v,l/w we obtain 
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Here the value wi 2 is defined differently depending on the sign of the product w\w 2 . For w\w 2 > we have 
w\ 2 = Wi + w 2 while in the opposite case wiw 2 < we define 1012 = max[|«;i|, \w 2 \}. We observe that in 
both cases the function F depends linearly on frequency and vanishes as soon as \w\ 2 \ exceeds \v\/2. 

At arbitrary values of v, wi and w 2 the integral (29) can be evaluated numerically. The corresponding 
result for the function F(v,W\,w 2 ) is presented in Fig. 2. 

Let us emphasize again that in the course of our analysis of the correlator ^(wi,^) we have ignored 
electron-electron interactions and also assumed that the external impedance Rs equals to zero, i.e. the 
voltage source is directly attached to the scatterer. In realistic experimental setups the latter condition 
may be violated, in which case the third cumulant can be substantially modified 21 by current fluctuations 
in the external environment (leads) . Within our formalism this effect can be accounted for by the external 
impedance Zs(oj). 
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Figure 2. Frequency dependence for the third cumulant of the current operator Sz{u\,u)2) evaluated by means of 
Eqs. (28)-(30) at T = for a short diffusive conductor with j3 = 1/3 and 7 = 2/15. 

4. CURRENT NOISE IN THE PRESENCE OF INTERACTIONS 

Let us now turn to the second cumulant (23) and explicitly account both for interactions and the external 
leads Rs- It is clear that fluctuations both in the scatterer and in the external resistor will affect the 
correlator (23), i.e. the latter should in general depend on both R and Rs- Let us - in addition to (23) 
- introduce another correlator S(t, t') defined by the same Eq. (23) in which one should substitute the 
current operator across the scatterer I — > I sc . The two correlators S(t,t') and S(t,t') are not independent. 
With the aid of the current conservation condition and performing the Fourier transformation with respect 
to t — if , one easily obtains the relation between these correlation functions: 



= tk(l+u 2 R%C 2 )S„- §(l+w 2 ^C 2 )wcoth^, (36) 
Hq re Zl 



where Rq = RRs/{R + Rs)- The second term is due the noise produced by the external resistor Rs which 
has to be subtracted in order to arrive at 5^. 

In general also the correlator <S W depends on both R and Rs- However, in the limit Rs 3> R the 
dependence on the shunt resistance is weak and can be neglected. In this case the interaction correction 
to the current noise spectrum depends only on the properties of the scatterer. Below we will present our 
results only in this limit. More general expressions can be found in Ref. 9. 

As before, we will assume that the condition g + gs 3> 1 is fulfilled throughout our calculation. Combin- 
ing the result (23) with the definition (24) and having in mind the relation (36) one can directly evaluate 
the correlator SS^, in the presence of interactions. In the important limit u>, eV, T -C l/RC one finds 

„ , \ 2(1-/3) , u „ ~ /UJ T \ , u + eV „ ~ /uj \ , w — eV 

S 2 (u) - 1 - ^ ^coth— + < 9ej(- + yj C oth-^+ < 9ej(--^jcoth-^- 

r°° r°° (l~c~ x / RC )(coseVy-coseV^on fil ,n l 

+ 7T7T e / ax ay — — — — — — — — 

Jo Sinn nTx smh nTy 

1 ^ (37) 



where 



sinli7rT(.T — y) sinh7rT(a; + y) 

I (V) = Y 5~- dt , (l- e -V^ )s ineW. (38) 

R 7T J smh nTt 



The result (37) covers several important regimes. In the limit of weak tunneling j3 — ► 1 and 7 — > wc find 
S 2{w) = e j(^ + V ) m th^^ + el(^-v)coth^- (39) 



e J 2T Ve / 2T 

Here the function I(V) is just the i" — V curve of a tunnel barrier modified by weak Coulomb blockade 
corrections. Eq. (39) can be viewed as a generalization of the result. 22 This formula also agrees with the 
result 23 derived in the limit g <C 1. 

It is worth pointing out that Eq. (37) can be reduced to the form (39) only in the case of tunnel barri- 
ers. Otherwise the interaction correction to the current noise cannot be obtained from the corresponding 
correction to the I — V curve except for equilibrium situations where the fluctuation-dissipation theorem 
(FDT) can be employed. In order to illustrate this point and to present our results beyond the tunneling 
limit let us split the function S into two parts: 

S(t,t') =S ni (t,t')+6S(t,t'), (40) 

where <S nl is the noninteracting contribution to the current noise 1 and SS is the correction due to electron- 
electron interactions inside the scatterer. Let us define the average voltage across the scatterer V = 
V X R/(R + Rs) and consider first the limit of relatively small voltages. At sufficiently large temperatures 
and/or frequencies wc find 

SS U = -^-, ]£T»gE c ,\eV\,\w\, (41) 



3i? 
3E 
R 

At lower temperatures and frequencies we obtain 



SS U = if M » T,gE c , \eV\. (42) 



5S u = -^\n^, if M,|eV|«T«0£ Ol (43) 
K q 1 



sSu = _2^j\ }ji gEc if T, \eV\ <C \u\ <C gEc- (44) 

Rq \U\ 

These results apply as long as either temperature or frequency exceeds the parameter gEc exp(— g/2). 

As we have already pointed out, the above expressions could also be anticipated from FDT combined 
with the results. 8 Indeed, in the limit of low voltages the current noise is described by the standard 
Nyquist formula. Hence, in order to satisfy FDT one should simply substitute the effective conductance 8 
into this formula. In this way one immediately reproduces Eqs. (41) and (43). 

Now let us turn to the case of relatively large voltages V where the shot noise becomes important. In 
this case the correction to the noise power spectrum is found to be proportional to the parameter 

^_ 27= E„T„(l-T) ( l-2T n)j (45) 



We obtain 



- 2(^-27)^1 , gE c 
R q \eV\ 

6S u = - ^~ 2 ^ )Ec , i£\eV\Z>T,gE c ,H (46) 



We note that this correction can be either negative or positive depending on the relation between the 
parameters /3 and 7. Thus, in contrast to the limit of low voltages (Nyquist noise), one cannot conclude 
that shot noise is always reduced by interactions. This reduction occurs only for conductors with relatively 
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Figure 3. Frequency dependence of the noise correlator 1S2 at T = and for various parameters of the system. 

low transmissions j3 > 27, while for systems with higher transmissions the net effect of the electron-electron 
interaction enhances the shot noise. In the important case of diffusive conductors one has (3 = 1/3, 7 = 2/15 
and, hence, (3 — 2-f — 1/15. In this case the shot noise is reduced by interactions. 

The above results have a transparent physical interpretation. At low voltages the power spectrum of the 
Nyquist noise is proportional to the system conductance oc J2 n Tn- Since in the presence of interactions the 
conductance acquires a correction proportional to f3, the interaction correction to the Nyquist noise should 
scale with the same parameter. On the other hand, shot noise is determined by the combination 5^_(T n — 
T%). Accordingly, the interaction correction to the shot noise power should consist of two contributions. 
One of them comes from T n and is again proportional to (3. Another contribution originates from the 
interaction correction to ^2 n T% which turns out to scale as 27. Since these two corrections enter with 
opposite signs we immediately arrive at the combination (45). 

We also point out that the third cumulant of the current operator for noninteracting electrons is known 5 
to be proportional to the parameters (3 and (3 — 27 respectively at low and high voltages. Following the 
same arguments as above one can anticipate that the interaction correction to the third cumulant should 
scale as [3 — 27 at low voltages and as (3 — 67 + 65 at high voltages, where 6 — J2 n ^nO- " ^»)/ S n ^ n 
Technically, in order to determine the interaction correction to the third cumulant 1S3 it is necessary to 
expand the exact effective action So (7) up to the fourth order in ip~ . Along with this interaction correction 
such expansion allows to determine the full frequency dependence of the fourth cumulant S4 in the absence 
of interactions. This rule also applies to higher cumulants, i.e. the lowest order interaction correction to 
the n-th current cumulant S n is determined by <S n +i for all values of n. This fact can also be proven by 
means of the renormalization group analysis recently developed in Refs. 11 and 12. 

5. CONCLUDING REMARKS 

In this paper we have presented a general path integral approach which allows to describe statistics of current 
fluctuations in mesoscopic coherent conductors at arbitrary frequencies and in the presence of interactions. 
This approach enables one to establish a complete expression for the effective action of coherent conductors 
described by an arbitrary - though energy independent - scattering matrix and to elucidate a profound 
relation between full counting statistics and electron-electron interaction effects in coherent mesoscopic 



conductors. Further extention of our technique to more complicated structures described by the energy 
dependent scattering matrices is possible and was recently worked out. 24 

Restricting ourselves to the non-interacting case, we have analyzed frequency dispersion of the third 
cumulant of the current operator. This dispersion was found negligible only in the case of tunnel junctions, 
while in a general case it turns out to be important in the frequency range comparable to eV. For instance, 
in the important case of diffusive conductors and at T = the quantity S3 changes by the factor 5 depending 
on whether relevant frequencies are below or above eV. For conductors with (3 < 27 even the sign of S3 
differs in these two limits. 

We have also analyzed the behavior of the second cumulant S2 in the presence of electron-electron 
interactions. We have demonstrated that Coulomb interaction decreases the Nyquist noise as one could 
anticipate already from the results 8 combined with FDT. In the case of the shot noise, the effect of electron- 
electron interactions turns out to be more complicated. The corresponding interaction correction was found 
negative for conductors with f3 > 27 and positive otherwise. 

Finally, let us point out that our predictions can be experimentally tested in various types of coherent 
mesoscopic conductors, such as, e.g., break junctions, quantum point contacts or short diffusive metallic 
bridges. 25,26 In all these systems both td and trc can be small enough in order to satisfy all the 
assumptions adopted here. For instance, in diffusive samples 25 ' 26 one finds 1/td of order few Kelvins and, 
hence, the condition eV < 1/tq is obeyed in a wide range of voltages V < 0.1 -j- 0.5 mV. We also note that 
the work 26 reports the experimental analysis of the frequency dispersion of shot noise 1 which is important 
in the same frequency range as that of the third cumulant studied here. 

This work is part of the Kompetenznetz "Funktionelle Nanostructuren" supported by the Landestiftung 
Baden- Wiirttemberg gGmbH and of the STRcP "Ultra-ID" supported by the EU. 
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